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Abstract 

The dijet cross section in photoproduction has been measured with the ZEUS 
detector at HERA using an integrated luminosity of 38.6 pb“^. The events were 
required to have a virtuality of the incoming photon, Q^, of less than 1 GeV^ 
and a photon-proton centre-of-mass energy in the range 134 < W^p < 277 GeV. 
Each event contains at least two jets satisfying transverse-energy requirements 
of > 14 GeV and > 11 GeV and pseudorapidity requirements of 
— 1 < ^ 2.4. The measurements are compared to next-to-leading-order 

QGD predictions. The data show particular sensitivity to the density of par- 
tons in the photon, allowing the validity of the current parameterisations to be 
tested. 
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1 Introduction 


In photoproduction at HERA, a quasi-real photon, emitted from the incoming positron, 
collides with a parton from the incoming proton. The photoproduction of jets can be 
classified into two types of process in leading-order (LO) Quantum Chromodynamics 
(QCD). In direct processes (Fig. 0(a)), the photon participates in the hard scatter via 
either boson-gluon fusion or QCD Compton scattering. The second class, resolved pro¬ 
cesses (Fig. |l|(b)), involve the photon acting as a source of quarks and gluons, with only 
a fraction of its momentum participating in the hard scatter. Measurements of jet cross 
sections in photoproduction [1-8] are thus sensitive to the structure of the photon and the 
proton, and to the dynamics of the hard sub-processes as calculated in perturbative QCD 
(pQCD). These jet cross sections can therefore be used in global hts to data to determine 
the parton densities in both the photon and proton. 

In the kinematic range of the measurements presented in this paper, the value of Xp, the 
fractional momentum at which partons inside the proton are probed, lies predominantly in 
the region between 0.01 and 0 . 1 . At these Xp values, the parton densities in the proton are 
constrained by measurements of the structure function, Fl’, [9] in deep inelastic scattering 
(DIS). The present measurements are directly sensitive to both the quark and the gluon 
content of the photon. The fractional momentum of the photon carried by the interacting 
parton, x-^, lies between 0.1 and 1. For x^ values above 0.5, the quark densities in 
the photon are not well constrained by data obtained from 77 * scattering in e’''e“ 
experiments [10-16]. The gluon density, to which jet photoproduction is directly sensitive 
at LO, is also poorly constrained by the F 2 data for all x^. The most recent measurements 
of F 2 from LEP extend up to an average scale of ~ 25 GeV. These, and higher scales, 
can be studied in jet production at HERA. 

The aim of the present investigation is to provide constraints, from data on dijet pho¬ 
toproduction, on the parton densities in the photon in the range 0.1 < < 1 and to 

probe the dynamics of the hard sub-processes. For this purpose, the dijet cross section 
is measured at high transverse energies where next-to-leading-order (NLO) QCD calcula¬ 
tions are expected to describe the data. In this kinematic region, where the effects of soft 
physics are suppressed and the parton densities in the proton are well known, the data 
can be used to test the validity of the current parameterisations of the parton densities in 
the photon. At high x-^, where the effects of the uncertainties in the photon structure are 
small, these data also provide a consistency check of the gluon distribution in the proton 
extracted from deep inelastic scattering. 

This analysis builds on the improved understanding of jet photoproduction and of com¬ 
parisons to NLO QCD calculations gained in previous analyses [2,4]. With an increase of 
a factor of six in luminosity, an extension of the kinematic region and reduced systematic 
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uncertainties, the measurements in this papeiQ have greatly improved precision compared 
to the previous ZEUS dijet measurement [2]. 


2 Theoretical framework 

Within the framework of pQCD, the dijet photon-proton cross section, da^p, can be 
written as a convolution of the proton parton density functions (PDFs), fp, and photon 
PDFs, fry, with the partonic hard cross section, d&ab- 


da^p = '^ dxpdx^fp{xp, fi‘^)f^{x^, fi^)ddab{xp, x^, jj,^) ■ (1 4ad)- 


For photoproduction cross sections measured in lepton-proton scattering, there is an ad¬ 
ditional convolution with the distribution of photons from the lepton beam. In the case 
of the direct cross section, the photon structure is replaced by a delta function aX x^ = 1. 
The scale of the process, /i, represents both the renormalisation scale, and factori¬ 
sation scale. Up, which are set equal for this study. The hadronisation correction, (5had, 
accounts for non-perturbative effects in the hnal state and can be estimated using Monte 
Carlo (MC) models for the parton cascade and fragmentation; it is, in general, a function 
of the variable being measured (see Section ||). 

The distribution of the dijet angle, 9*, in the parton-parton centre-of-mass frame is directly 
sensitive to the form of the matrix elements, and hence to the partonic hard cross section. 
For massless partons, the centre-of-mass scattering angle is given by 


cos 9* 


tanh 





where and are the pseudorapidities in the laboratory frame of the two jets of 
highest transverse energy. Only the absolute value of cos 9* can be determined because the 
originating parton cannot be identihed. The variable cos 9* is invariant under the different 
boosts along the beam axis arising from the spectrum of incoming parton momenta. This 
minimises the sensitivity of the differential cross section, da/dcos9*, to the momentum 
density distribution of the partons in the photon and proton. 

^ After submission of this paper, a similar study from the HI collaboration became available [17]. 
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For jets of transverse energy of more than 6 GeV, it has been shown [5] that samples 
of events enriched in either direct or resolved photon processes have very different an¬ 
gular distributions. The cross section for the sample enriched in resolved photon events 
increases more rapidly at high cos 6* than that in direct photon events. This is expected 
at both LO and NLO QCD [18]; both predictions give a good description of the data. 
The different angular dependence of the cross sections can be explained in terms of the 
dominant propagators in the respective samples. In direct events, the dominant processes 
(mostly boson-gluon fusion) have a spin-| quark propagator and the angular dependence 
is approxiamtely oc (1 — |cos6**|)“^. In resolved events (e.g. qg —> qg and gg —> gg), the 
dominant processes have a spin-1 gluon propagator, which has an angular dependence 
oc (1 — |cos6**|)“^. 

To probe the structure of the photon, it is desirable to measure cross sections as functions 
of variables that are sensitive to the spectrum of incoming parton momenta, such as the 
pseudorapidity of the jets or the momentum fraction, x^. Since is not directly measur¬ 
able, an observable, x'^^, which is the fraction of the photon’s momentum participating 
in the production of the two highest transverse-energy jets (and is equal to x^ for partons 
in LO QCD), is introduced [6]: 


X 


obs 

7 


‘^yEe 


Here and are the transverse energies of the two jets in the laboratory frame 
and y is the fraction of the positron’s energy, T’g, carried by the photon in 
the proton rest frame. The quantity x^^ is a particularly useful variable with which to 
discriminate between various photon PDFs. 


3 Experimental conditions 

The data were collected during the 1996 and 1997 running periods, when HERA operated 
with protons of energy Ep = 820 GeV and positrons of energy E^ = 27.5 GeV, and 
correspond to an integrated luminosity of 38.6 ± 0.6 pb~^. A detailed description of the 
ZEUS detector can be found elsewhere [19,20]. A brief outline of the components that 
are most relevant for this analysis is given below. 

The high-resolution uranium-scintillator calorimeter (CAL) [21] consists of three parts: 
the forward, the barrel and the rear calorimeters. Each part is subdivided transversely 
into towers and longitudinally into one electromagnetic and either one (in the rear) or two 
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(in the barrel and forward) hadronic sections. The smallest subdivision of the calorime¬ 
ter is called a cell. The CAL relative energy resolutions, as measured under test-beam 
conditions, are 0.18/V^ for electrons and 0.35/v^ for hadrons {E in GeV). 

Charged particles are measured in the central tracking detector (CTD) [22], which op¬ 
erates in a magnetic held of 1.43 T provided by a thin superconducting coil. The CTD 
covers the polar-angle^ region 15° < 9 < 164°. The relative transverse-momentum resolu¬ 
tion for full-length tracks is o-p^/pr = 0.0058pT © 0.0065 © 0.0014 /pt with in GeV. 
Tracking information along with energy deposits in the CAL were used to measure the 
transverse energy and direction of jets as described in detail in Section |^. 

The luminosity was measured from the rate of the bremsstrahlung process e'^p —>■ e'^'yp, 
where the photon was measured in a lead-scintillator calorimeter [23] placed in the HERA 
tunnel at Z = —107 m. 

A three-level trigger system was used to select events online [4,20]. At the third level, 
a cone algorithm was applied to the CAL cells and jets were reconstructed using the 
energies and positions of these cells. Events with at least two jets, each of which satished 
the requirements that the transverse energy exceeded 4 GeV and the pseudorapidity was 
less than 2.5, were accepted. 


4 Definition of the cross section 


The kinematic region for this study is the photoproduction region, dehned as < IGeV^, 
with a photon-proton centre-of-mass energy, HT^p, in the range 134 GeV to 277 GeV. Each 
event is required to have at least two jets reconstructed with the fcr cluster algorithm [24] 
in its longitudinally invariant inclusive mode [25], with at least one jet having transverse 
energy greater than 14 GeV and another greater than 11 GeV. The jets are required to 
satisfy — 1 < < 2.4, an extension of the pseudorapidity range by 0.4 units in the 

forward direction over the previous analysis [2], thereby increasing the sensitivity of the 
measurement to resolved photon processes. 

Cross sections are presented as a function of and The cross sections for 

jet variables are symmetrised [26] with respect to the pseudorapidities of the two jets. 
Each event, therefore, contributes twice to the cross section. The cross sections have been 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the “forward direction”, and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity 
is defined as ry = — ln(tan |), where the polar angle, 6, is measured with respect to the proton beam 
direction. 
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determined for regions enriched in direct and resolved photon processes by requiring 
to be greater than 0.75 or less than 0.75, respectively. 

Additional kinematic constraints were applied to the measurement of the cross section as 
a function of |cos0*| to remove biases imposed by the other requirements. For a given 
centre-of-mass energy, events at high |cos0*| have small scattering angles and thus lower 
To study the | cos6**I distribution up to |cos0*| = 0.8 without bias from the 77^* 
requirements, a cut on the dijet mass of > 42 GeV was applied. The dijet mass is 
dehned in terms of the two jets with highest transverse energy as 


Afjj = sj [cosh(r7Fti - _ cos ((/>)®*2 _ 0 jet 2 )] ^ 

where and are the azimuthal angles of the two jets. For jets back-to-back in 
azimuthal angle and of equal transverse energy, 77;^*, the dijet mass is related to the 
scattering angle by 


MJj 


2E^^ 

■\/l — |cos6'*P 


( 1 ) 


When the minimum 77;^* is taken to be 12.5 GeV, the average of the minimum transverse 
energies of the two jets, the requirement on the minimum dijet mass up to a given value 
of the scattering angle can be deduced from Eq. (|T[). Simulation studies show that the 
choice of cut dictated by this approximation does indeed eliminate any bias from the 
choice of transverse-energy cuts [27]. A further cut on the boost of the dijet system in the 
laboratory frame, -|-r7j®*^)/2, of 0.1 < f] < 1.3 was also applied. This ensures that 

the phase space is uniform as a function of |cos6**|, so that any shape seen in the measured 
distributions is attributable to the dynamics and not biased by the cuts imposed. 


5 Monte Carlo models 

The acceptance and the effects of detector response were determined using samples of 
simulated events. The programs Herwig 6.1 [28] and Pythia 6.1 [29], which implement 
the leading-order matrix elements, followed by parton showers and hadronisation, were 
used. The Herwig and Pythia generators differ in the details of the implementation 
of the leading-logarithmic parton-shower models. They also use different hadronisation 
models: Herwig uses the cluster [30] model and Pythia uses the Lund string [31] model. 
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Direct and resolved events were generated separately. For all generated events, the ZEUS 
detector response was simulated in detail using a program based on Geant 3.13 [32]. 

Parameters tuned to HERA data [33] were used in the generation of the Herwig sample. 
The GRV-LO [34] and CTEQ4L [35] set of PDFs were used for the photon and proton, 
respectively. For the Pythia generator, the parameters were chosen to be consistent with 
hts to jet data from both HERA and LEP [36]. Here, the SaS-2D [37] and GRV94-LO [38] 
set of PDFs were used for the photon and proton, respectively. Multiparton interactions 
(MPI) were also included with a minimum transverse momentum of the secondary scatter 
of 2.0 GeV [39]. However, at the high transverse energies studied here, the effects of the 
“underlying event” are small: models with or without MPI describe the data equally well. 


6 NLO QCD calculations 


Many calculations of jet photoproduction at NLO exist [40-46], all of which have been 
compared with each other and agree to within (5 — 10)% [46,47]. The calculation used 
here is that of Frixione and Ridolh [40,41], which employs the subtraction method [48] 
for dealing with the collinear and infra-red divergences. The number of flavours was set 
to 5 and the renormalisation and factorisation scales, /i, were set to half the sum of the 
transverse energies of the hnal-state partons, Et/2. Two different parameterisations of 
the photon parton density were used|^: GRV-HO [34] and AFG-HO [51]. The parton den¬ 
sities in the proton were parameterised using CTEQ5M1 [35]; the value as{Mz) = 0.118 
used therein was adopted for the central prediction. The parameterisation MRST99 
{as{Mz) = 0.1175) was also considered. Parameter settings in the NLO calculation were 


varied to test the stability of the theoretical predictions, as discussed in Section 10 


The NLO QCD predictions were corrected for hadronisation effects using a bin-by-bin 
procedure according to da = da^^^ ■ where is the cross section for par- 

tons in the hnal state of the NLO calculation. The hadronisation correction factor, 
Chad = (1 + dhad)~^ was dchued as the ratio of the dijet cross sections before and af¬ 
ter the hadronisation process, Gy,ad = The value of Chad was taken as 

the mean of the ratios obtained using the Herwig and Pythia predictions. The hadro¬ 
nisation correction, dh ad , was generally below 10% in each bin except at the edges of phase 
space (see Tables M)- There are signihcant migrations in for >0.7; however, 
the migrations to lower values are small. 


^ The GS96-HO [49] parameterisation was not used because it has been shown that the available code 
does not reproduce the published results [50]. 
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7 Energy corrections 


Kinematic variables and jets were reconstrncted using a combination of track and calorime¬ 
ter information that optimises the resolution of reconstructed kinematic variables [52], 
The selected tracks and calorimeter clusters are referred to as Energy Flow Objects 
(EFOs). 

The addition of track information to the CAL information reduces the sensitivity to en¬ 
ergy losses in inactive material in front of the CAL, and exploits the good momentum and 
angular resolution of the tracking for low-momentum particles. The energies of particles 
for which no track information was available (e.g. neutral particles), or for which the 
calorimeter energy resolution was better than that of the tracking (e.g. at the highest 
energies), were measured using CAL information. These energies were corrected for losses 
in the inactive material as discussed in a previous publication [2]. Conservation of en¬ 
ergy and momentum in neutral current (NC) DIS events was exploited to determine the 
required energy corrections [53] by balancing the scattered positron with the hadronic 
hnal state. This was performed independently for data and simulated event samples. 
The EFOs thus corrected were used both to reconstruct jets and to determine kinematic 
variables. Comparisons of kinematic variables for data and simulated events led to the 
assignment of a 1% correlated systematic uncertainty in the transverse jet energies and in 
the hadronic variables [53]. The improved precision in this uncertainty compared with the 
previous measurement [2] was obtained from an increased data sample, better selection 
requirements and improved parameterisations of the energy losses. In the overlapping 
kinematic region, the total cross section measured here is 6% lower than the previous 
measurement [2], but within the quoted uncertainties arising from the uncertainty in the 
jet energy scale. 


8 Event selection 

After applying the energy corrections described in Section |^, dijet events were selected 
offline by using the following procedures and cuts designed to remove sources of back¬ 
ground: 

• the kx clustering algorithm was applied to the corrected EFOs. Events were selected 

in which at least two jets were found with —1 < < 2.4, > 14 GeV and 

> 11 GeV; 

• to remove background due to proton beam-gas interactions and cosmic-ray showers, 
the longitudinal position of the reconstructed vertex was required to be in the range 
I ^vertex | ^ 40 Cm, 
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• to remove background due to charged current DIS events and cosmic-ray showers, a 
cut on the relative transverse momentum of pr/y/Er < l.SVGeV was made, where px 
and Et are, respectively, the measured transverse momentum and transverse energy 
of the event; 

• NC DIS events with a scattered positron candidate in the CAL were removed by 
cutting [7] on the inelasticity, y, which is estimated from |/e = 1 — ^(1 — cos6*'), 
where E'^ and 9'^ are the energy and polar angle, respectively, of the scattered positron 
candidate. Events were rejected if pe < 0.85; 

• the reqnirement 0.2 < yjB < 0.85 was imposed, where yjB is the estimator of y 
measured from the CAL energy deposits according to the Jacquet-Blondel method [54]. 
The upper cut removed NC DIS events where the lepton was not identihed and which 
therefore have a value of |/jb close to 1. The lower cut removed proton beam-gas events 
which have a low value of ?/jb. These requirements on ?/jb restrict the photon-proton 
centre-of-mass energy to be in the range 134 to 277 CeV. 

The cuts on y^. and yjB reduced the background from DIS events to less than 0.5% and 
restricted the range of the virtuality of the exchanged photon to < 1 CeV^, with a 
median value of about 10“^ CeV^. After these requirements, 62573 events with two or 
more jets remained in the sample; 2919 of these events had a third jet of transverse energy 
greater than 11 CeV in the range —1 < < 2.4. For the measurements of the cross 

section as a function of cos0*, the following additional requirements were imposed: 

• the dijet scattering angle was restricted to be in the range jcos^*! < 0.8; 

• the invariant mass of the two jets of highest transverse energy was required to satisfy 
Afjj > 42 CeV; 

• the boost of the dijet system was required to be in the range 0.1 < fj < 1.3. 

These cuts reduced the sample to 10811 events. 


9 Event characteristics 

For the high transverse energies studied here, it has been shown previously [2,4] that the 
transverse energy flow around jets is generally well described by the simulation. This 
agreement is maintained in the more-forward region studied here, with Pythia giving a 
description similar to that of Herwig. 

Kinematic distributions in the data are compared to the two simulation programs in 
Fig. i The simulated distributions were ht to the data distribution in shown in 
Fig. 0(a), by varying the fractions of direct and resolved processes and minimising the 



The simulations generally describe the data well for all variables, although some 
discrepancies are seen. The shape of the distribution in yj-Q is better described by the 
Herwig prediction, although the description by the Pythia simulation is adequate. Both 
simulations have similar distributions for the transverse energies of the jets and describe 
the data well. However, neither gives a good description of the pseudorapidity of the jet of 
highest transverse energy, whereas the pseudorapidity distribution of the second jet is well 
described by Herwig. Figure also shows that according to the Herwig simulation, the 
proportion of resolved photon events decreases with increasing and and increases 
with increasing |/jb and 

Since Herwig gives a better overall description of the data than Pythia, it was chosen 
as the primary MC generator to correct the data. The correction was performed using the 
bin-by-bin method, in which the correction factor, as a function of an observable (P in a 
given bin i, is Ci{0) = The variable is the number of events 

passing the kinematic requirements on the hadronic hnal state described in Section ^ 
and is the number of reconstructed events passing the selection requirements as 

detailed in Section Both numbers were computed using the MC generators described 
in Section For distributions as a function of the correction factors lie in the range 
0.85-1.25. 


10 Experimental and theoretical uncertainties 

10.1 Experimental uncertainties 

The results of a detailed analysis of the possible sources of systematic uncertainty are listed 
below. Typical values for the systematic uncertainty are quoted for the cross sections as 
a function of 

• varying the measured jet energies by ±1% in only the simulated sample, in accordance 
with the uncertainty in the jet energy scale, gave an uncertainty of =f 5%, increasing 
with 

• correcting the data with Pythia instead of the Herwig generator gave an uncertainty 
within ±9% and typically ±4%; 

• changing the cuts on and in both the data and simulated samples by the 
value of the average resolution (~ 9%) gave an uncertainty of ±5%; 

• changing the cuts on in both the data and simulated samples by the value of the 
resolution (~ 0.03 at low yjB and ~ 0.05 at high |/jb) gave an uncertainty of less than 
± 2 %; 
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• changing the cuts on both the data and simulated samples by the value of 

the resolution (~ 0.04) gave an uncertainty of ±0.5%; 

• varying the cuts to remove DIS and beam-gas backgrounds in both data and simulated 
samples gave a total uncertainty of less than ±1%. 

In addition to the above, the cuts made to evaluate the cross section as a function of 
I cos 0*1 also lead to sources of systematic uncertainty. The following were evaluated, with 
typical uncertainties quoted: 

• changing the cuts on Tfjj in both the data and simulated samples by the value of the 
average resolution (~ 8%) gave an uncertainty of ±5%; 

• changing the cuts on fj in both the data and simulated samples by the value of the 
resolution (~ 0.04) gave an uncertainty of ±2%; 

• changing the cuts on in both the data and simulated samples by the value of 
the resolution (~ 0.04) gave an uncertainty of ±4%. This systematic uncertainty also 
contributes to the measurements of the cross-sections da/dE}^^^ and da/drf'^^'^. 

The uncertainty in the cross sections due to the jet energy scale uncertainty is correlated 
between bins and is therefore displayed separately as a shaded band in the hgures. All 
other systematic uncertainties were added in quadrature. In addition, an overall normal¬ 
isation uncertainty of 1.6% from the luminosity determination is not included in either 
the hgures or tables. 

As a cross check, cross sections were obtained with an iterative matrix-unfolding technique 
[55], using Bayes’ theorem. The resultant cross sections were found to be consistent with 
those using the standard bin-by-bin procedure [27]. 

10.2 Theoretical uncertainties 

The NLO QCD predictions for the dijet cross section are affected by the systematic 
uncertainties listed below. Typical values for the systematic uncertainty are quoted for 
the cross sections as a function of 

• the uncertainty due to terms beyond NLO, estimated by varying /i between Ex and 
A'r/4, is ±(10 — 20)%, increasing with decreasing 

• the uncertainty due to the hadronisation correction, estimated as half the spread 
between the (5had values obtained using the Herwig and Pythia models, is ±(2 —3)%; 

• the uncertainty due to the value of as{Mz), estimated by repeating the calculations 
using the CTEQ4 series of PDFs determined using values of as{Mz) = 0.113, 0.116 
and 0.119, is ±(5 - 8)%. 
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The above systematic uncertainties are largely independent of the choice of photon PDF 
and were added in quadrature to give the total uncertainty on the predictions in each 
case. Differences between parameterisations of the photon and proton parton densities 


are discussed in the comparison to the measnred data in Section 11 


11 Results 

11.1 Probing the matrix elements 

The dijet cross section as a function of jcos^*!, is given in Table and shown in Fig. |^. 
The data are shown separately for < 0.75 (Fig. |(a)) and > 0.75 (Fig. |(b)) 
and compared to NLO predictions. For < 0.75, the measured cross section lies above 
the NLO prediction using GRV-HO for the photon PDF by an average of (10 — 15)%. 
Considering the theoretical and experimental uncertainties, both of (5 — 10)%, the NLO 
prediction gives a reasonable description of the data. The predictions using the AFG-HO 
parameterisation for the photon give a lower cross section than that of GRV-HO, and are 
thus around (20 — 25)% lower than the data. For > 0.75, the NLO prediction is in 
agreement with the measured cross section. 

In Fig. |^(c), the shapes of the data and NLO distribntions are compared. The predictions 
give a generally good description of the data; the shapes of the predictions when using 
the GRV-HO and AFG-HO parameterisations are similar. The data at low rise more 
rapidly at high |cos0*| than those at high This is consistent with a difference in the 
dominant propagators, as observed in a previous publication [5]; this is seen here at higher 
energies and masses. The agreement in shape of these distributions, which are sensitive 
to the matrix elements, demonstrates that also in this high-mass region the dynamics of 
the short-distance process is understood. 

11.2 Cross sections in and 

Measurements of the dijet cross section as a function of are given in Tables and 
^ and shown in Figs. ^-|^, in different regions of the pseudorapidities of the two jets for 
> 0.75 and < 0.75. The ratios of the cross sections to theory are shown in Figs. |] 
and 1^. In Fig. the measurement for > 0.75 extends to transverse energies of ~ 70 
GeV, extending the region measured previously [2]. In general, the overall description of 
the data by the predictions is reasonable. In particular, when the jets are produced in 
the region 1 < 77 ^®*^ <2.4 and 0 < < 1, the cross section in Fig. ^ falls three orders of 

magnitude and is well described by the NLO calculation. When both jets are produced 
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in the region 1 < < 2.4, the NLO prediction lies below the data at low transverse 

energy, althongh both the experimental and theoretical nncertainties are sizeable. In 
this region, the hadronisation corrections are signihcant bnt do not acconnt for all the 
differences. This region is near the edge of phase space dne to the cnts applied. 

At low (Figs. | and 0 ). the data are also generally well described by the NLO predic¬ 
tions, althongh a difference in shape is seen, with the predictions lying above the data at 
low transverse energy and below for > 20 GeV. The predictions nsing AFG-HO are 
nniformly abont 15% below those of GRV-HO for the entire range of transverse energies. 

The psendorapidities of the two jets are sensitive to the momentnm distribntions of the 
incoming partons. The cross section is measnred as a fnnction of the psendorapidity of 
one of the jets, in different regions of psendorapidity of the other, and is shown at high 
^obs (^pig_ l^j and low (Fig. The NLO predictions give a good description of the 
data except for — 1 < < 0 for low where the data are at or below the lower 

edge of the scale-nncertainty band. The predictions nsing AFG-HO lie abont (10 — 15)% 
below those of GRV-HO. 


11.3 Testing the current parameterisations of the photon PDF 


The cross sections and ratios of data and theory as a fnnction of in regions of increasing 


transverse energy are shown in Figs. pI)|-[T^. The predictions lie signihcantly above the 
data nsing the GRV-HO parameterisation at the lowest valnes (14 < < 17 GeV) 

of transverse energy for > 0.5, bnt are increasingly below the data for valnes larger 
than 17 GeV. This trend with transverse energy is stronger for < 0.8, as can be seen 
in Fig. 1^. Given the uncertainties, the data and predictions are consistent except in the 
region of lowest transverse energy for > 0.5. The dominant theoretical uncertainty, 
estimated from the variation of the renormalisation and factorisation scales, arises from 
the higher-order contributions not present in an NLO calculation. The inclusion of higher- 
order contributions would have to result in a signihcant change of shape of the distribution 
as a function of both the transverse energy and if it were to describe the data. The 
data in Fig. |T^ are also compared to the NLO prediction using the AFG-HO photon 
parameterisation; the ratio of data to the theory is shown in Fig. The prediction 
agrees with the data in the region of lowest transverse energy, but is below the data for 
the higher Ei^^^ bins. The predictions using AFG-HO are similar in shape to those using 
GRV-HO but are (10 — 15)% lower. Figures and |T^ also show the predictions using the 
MRST99 PDFs in the proton. The differences between the predictions with GTEQ5M1 
and MRST99 are everywhere less than 5%. 
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11.4 Discussion 


To improve the understanding of the features of the cross section in different regions of 
transverse energy, the sensitivity of the above comparisons to the value of the cut on 
the second jet has been studied. Starting at a minimum of 11 GeV, the cut on the 
second jet was raised in both data and theory for the region 25 < < 35 GeV; the 

results are shown in Fig. |^. With increasing the data fall, as expected; the 

trend is well reproduced by the Herwig simulation, which includes leading-logarithmic 
parton showers. The prediction from Herwig is normalised to the data in the first bin 


in Fig. |T^(a). The LO prediction (not shown) for this cross section is flat, since only 
two partons are emitted, which must have equal transverse energies by conservation of 
energy. The predictions of the shape of this distribution from 0{aa‘l) QGD are therefore 
the lowest non-trivial order predictions. The predictions from NLO QGD are shown; they 
fall less rapidly at low more rapidly at high ^Ft 2 ,cut data. Frixione 

and Ridolh [41] have shown that, when the requirements on the minimum transverse 
energy of the two jets are similar, the NLO calculation is infrared sensitive. This has 
been investigated further by considering the cross section for regions of high and low 
For > 0.8, shown in Fig. |^(b), the data and NLO QGD converge for low 
both being reasonably insensitive to the cut and similar in shape. From Fig. |T^, it can 
be seen that for a cut on the first jet of 25 < < 35 GeV, the cut on the jet of 

lower transverse energy has to be below 21 GeV for the NLO predictions to agree with 
the data. Figure |^(c) shows the region < 0.8; the predictions lie below the data at 
low but within the theoretical uncertainties. The prediction using AFG-HO is 

about (10 — 15)% below that of GRV-HO, but is similar in shape and is therefore just 
compatible with the data. 


The difference in the behaviour of the data and the calculations in Fig. |T^ implies that 
there is a significant dependence on comparisons between the measurements 

in the previous section and NLO QGD. By adjusting ^;Ft 2 ,cut each 

range, it would be possible to achieve agreement between the NLO prediction and the 
data. However, this seems to be a somewhat arbitrary procedure. 


The agreement with theory at high and high transverse energy, where the dependence 
on the photon structure is small, demonstrates a consistency between these data and the 
gluon distribution in the proton extracted from DIS data. Further discrimination between 
the current PDFs is currently not possible given the large uncertainties in the theory at 
low transverse energies and both the experimental and theoretical uncertainties at higher 
transverse energies. However, the data shown here signihcantly constrain the parton 
densities in the photon. These constraints would be made more stringent with improved 
higher-order, or resummed, calculations. 
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12 Conclusions 


The dijet cross section in photoproduction has been measured in the kinematic region 
Q2 < 1 GeV^, 0.2 < y < 0.85, > 14 GeV, > 11 GeV and -1 < < 2.4. 

In the high-mass region dehned by > 42 GeV and 0.1 < fj < 1.3, the dijet angu¬ 
lar distribution of the data is well reproduced by the NLO predictions, indicating that 
the dynamics of the short-distance process is understood. Over the wider region, the 
measurements are compared with NLO predictions using different parameterisations for 
the parton densities of the photon. The data fall less steeply with increasing transverse 
energy than do the NLO QCD predictions, and show sensitivity to the parton densities 
of the photon. Neither the AFG-HO nor the GRV-HO parameterisation, convoluted with 
the NLO matrix elements, fully describes all features of the data. There is agreement 
with theory at high and high transverse energy, where the dependence on the photon 
structure is small, which represents a consistency check of the gluon distribution in the 
proton extracted from deep inelastic scattering. The data at low signihcantly con¬ 
strain the parton densities in the photon; future parameterisations of the photon PDFs 
should take them into account. These constraints would be made more stringent were 
improved higher-order or resummed calculations available. 
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•^7 

< 0.75 





0 . 0 , 

0.1 

44.5 

± 3.5 

+1.6 

-3.8 

+2.7 

-2.0 

0.994 

± 

0.010 

0 . 1 , 

0.2 

40.1 

± 3.3 

+3.2 

-2.4 

+2.3 

-1.3 

1.027 

± 

0.021 

0 . 2 , 

0.3 

48.1 

± 3.6 

+2.5 

-3.3 

+2.4 

-2.1 

1.006 

± 

0.001 

0 . 3 , 

0.4 

54.5 

± 3.7 

+4.4 

-0.9 

+2.0 

-2.5 

1.022 

± 

0.017 

0 . 4 , 

0.5 

83.9 

± 4.7 

+2.7 

-5.4 

+4.5 

-3.6 

1.008 

± 

0.014 

0 . 5 , 

0.6 

115.7 

± 5.5 

+1.4 

-4.1 

+4.8 

-4.5 

1.028 

± 

0.020 

0 . 6 , 

0.7 

205.5 

± 7.6 

+14.7 

-4.7 

+8.1 

-10.7 

1.024 

± 

0.010 

0 . 7 , 

0.8 

406.9 

± 11.5 

+9.5 

-32.7 

+16.3 

-15.3 

1.031 

± 

0.001 




™obs 

•^7 

> 0.75 





0 . 0 , 

0.1 

141.8 

± 7.1 

+4.0 

-8.1 

+5.5 

-5.2 

0.973 

± 

0.003 

0 . 1 , 

0.2 

142.1 

± 7.1 

+2.7 

-10.0 

+5.7 

-5.6 

0.969 

± 

0.010 

0 . 2 , 

0.3 

155.3 

± 7.3 

+7.9 

-6.8 

+6.5 

-6.4 

0.958 

± 

0.016 

0 . 3 , 

0.4 

182.7 

± 7.9 

+5.2 

-7.7 

+7.3 

-6.5 

0.970 

± 

0.007 

0 . 4 , 

0.5 

233.8 

± 8.8 

+9.6 

-10.7 

+8.8 

-8.9 

0.964 

± 

0.006 

0 . 5 , 

0.6 

283.9 

± 9.7 

+9.1 

-10.8 

+9.4 

-9.9 

0.963 

± 

0.008 

0 . 6 , 

0.7 

439.7 

± 12.2 

+9.2 

-17.2 

+16.3 

-17.8 

0.946 

± 

0.020 

0 . 7 , 

0.8 

686.6 

± 15.9 

+26.0 

-18.3 

+24.8 

-22.6 

0.933 

± 

0.020 


Table 1: Measured cross sections as a function of | cos0 *| for < 0.75 and 
> 0 . 75 . The statistical, systematic and jet energy scale, Aes , uncertainties are 
shown separately. The multiplicative hadronisation correction applied to the NLO 
prediction is shown in the last column. The uncertainty shown for the hadronisation 
correction is half the spread between the values obtained using the Herwig and 
Pythia models. 
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bin ( GeV ) 

da/dE^^^ 

^stat 

-^syst 

^ES ( pb / GeV ) 

Ghad 



-1 < 

^jetl,2 ^ Q 


14 , 17 

12.2 

± 0.6 

+ 3.6 

- 1.1 

+0.6 

-0.5 

0.794 ± 0.034 

17 , 21 

3.0 

± 0.3 

+ 0.4 

- 0.6 

+0.1 

-0.1 

0.757 ± 0.046 


0 < 

^jetl < _1 < 

^jet2 ^ Q 


14 , 17 

34.3 

± 0.7 

+ 2.2 

- 1.6 

+1.3 

-1.0 

0.909 ± 0.024 

17 , 21 

16.0 

± 0.4 

+ 0.5 

- 2.0 

+0.7 

-0.5 

0.907 ± 0.016 

21 , 25 

4.0 

± 0.2 

+ 0.1 

- 0.5 

+0.1 

-0.2 

0.863 ± 0.007 

25 , 29 

0.74 

± 0.08 

+0.08 

-0.17 

+0.02 

-0.03 

0.834 ± 0.013 



0 < 

^jetl,2 ^ ^ 


14 , 17 

51.4 

± 1.1 

+0.4 

-5.7 

+1.6 

-0.9 

0.954 ± 0.030 

17 , 21 

26.1 

± 0.7 

+0.7 

-2.1 

+0.9 

-1.0 

0.974 ± 0.017 

21 , 25 

12.3 

± 0.5 

+0.3 

-1.5 

+0.4 

-0.4 

0.973 ± 0.015 

25 , 29 

6.3 

± 0.3 

+0.1 

-0.8 

+0.3 

-0.2 

0.958 ± 0.021 

29 , 35 

1.7 

± 0.1 

+0.03 

-0.2 

+0.1 

-0.1 

0.953 ± 0.003 

35 , 41 

0.56 

± 0.08 

+0.02 

-0.08 

+0.03 

-0.02 

0.920 ± 0.006 

41 , 48 

0.122 

± 0.037 

+0.031 

-0.014 

+0.003 

-0.010 

0.872 ± 0.044 


1 < 

^jetl < 2 . 4 , -1 < 

^jet2 ^ Q 


14 , 17 

28.1 

± 0.6 

+3.4 

-0.7 

+1.0 

-1.1 

0.913 ± 0.015 

17 , 21 

15.8 

± 0.4 

+0.3 

-1.2 

+0.6 

-0.6 

0.922 ± 0.028 

21 , 25 

5.4 

± 0.2 

+0.2 

-0.1 

+0.2 

-0.1 

0.922 ± 0.012 

25 , 29 

1.6 

± 0.1 

+0.2 

-0.1 

+0.1 

-0.1 

0.923 ± 0.012 

29 , 35 

0.41 

± 0.05 

+0.02 

-0.05 

+0.01 

-0.02 

0.862 ± 0.014 

35 , 41 

0.043 

± 0.014 

+0.013 

-0.013 

+0.004 

-0.003 

0.827 ± 0.004 


Table 2: Measured cross section as a function of for events with > 0 . 75 . The 
measurement is divided into six regions of the pseudorapidities of the jets. For further 
details, see the caption to Table [J. 
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bin ( GeV ) 

da/dE}^^^ 

-^stat 

-^syst 

Aes ( pb / GeV ) 

Chad 


1 < 

< 2 . 4 , 0 < < 1 


14 , 17 

38.9 

± 0.7 

+ 0.4 

- 2.2 

+ 1.1 

- 1.1 

0.990 ± 0.010 

17 , 21 

23.1 

± 0.4 

+ 0.2 

- 1.2 

+ 0.8 

- 0.9 

0.978 ± 0.011 

21 , 25 

11.2 

± 0.3 

+ 0.1 

- 0.6 

+ 0.4 

- 0.3 

0.967 ± 0.017 

25 , 29 

5.3 

±0.2 

+ 0.03 

- 0.2 

+ 0.2 

- 0.2 

0.969 ± 0.015 

29 , 35 

2.1 

±0.1 

+ 0.01 

- 0.06 

+ 0.1 

- 0.1 

0.962 ± 0.022 

35 , 41 

0.87 

± 0.07 

+ 0.06 

- 0.05 

+ 0.06 

- 0.03 

0.971 ± 0.016 

41 , 48 

0.34 

± 0.04 

+ 0.01 

- 0.02 

+ 0.02 

- 0.02 

0.982 ± 0.009 

48 , 55 

0.13 

± 0.03 

+ 0.01 

- 0.03 

+ 0.01 

- 0.01 

0.954 ± 0.023 

55 , 65 

0.035 

±0.011 

+ 0.000 

- 0.014 

+ 0.002 

- 0.002 

0.974 ± 0.032 



1 < 7 ] 

etl ,2 ^ 2.4 


14 , 17 

1.5 

±0.2 

+ 0.4 

- 0.4 

+ 0.0 

- 0.1 

1.302 ± 0.225 

17 , 21 

4.0 

±0.2 

+ 0.8 

- 0.5 

+ 0.2 

- 0.2 

1.190 ± 0.032 

21 , 25 

5.0 

± 0.3 

+ 0.3 

- 0.1 

+ 0.2 

- 0.2 

1.083 ± 0.032 

25 , 29 

3.3 

±0.2 

+ 0.3 

- 0.2 

+ 0.1 

- 0.1 

1.017 ± 0.030 

29 , 35 

1.9 

±0.1 

+ 0.1 

- 0.1 

+0.0 

- 0.1 

1.010 ± 0.012 

35 , 41 

0.93 

±0.10 

+ 0.04 

- 0.03 

+ 0.08 

- 0.02 

0.992 ± 0.009 

41 , 48 

0.41 

± 0.06 

+ 0.06 

- 0.02 

+ 0.02 

- 0.03 

0.982 ± 0.018 

48 , 55 

0.23 

± 0.05 

+ 0.03 

- 0.04 

+ 0.01 

- 0.01 

1.000 ± 0.018 

55 , 65 

0.033 

± 0.015 

+ 0.029 

- 0.008 

+ 0.002 

- 0.002 

0.999 ± 0.004 

65 , 75 

0.041 

± 0.017 

+ 0.004 

- 0.010 

+ 0.002 

- 0.002 

0.987 ± 0.014 


Table 2: (cont.) 
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bin ( GeV ) 

da/dE^^^ 

^stat 

-^syst 

^ES ( pb / GeV ) 

Ghad 




-1 < 

^ jetl ,2 ^ Q 


14 , 17 

0.30 


± 0.08 

+ 0.07 

- 0.04 

+ 0.01 

- 0.01 

1.116 ± 0.130 



0 < 

^jetl < _1 < 

j ^ jet 2 ^ g 


14 , 17 

6.2 


± 0.3 

+ 0.9 

- 0.3 

+ 0.2 

- 0.1 

1.006 ± 0.003 

17 , 21 

2.1 


± 0.1 

+ 0.1 

- 0.4 

+ 0.1 

- 0.1 

1.014 ± 0.012 

21 , 25 

0.23 


± 0.04 

+ 0.05 

- 0.06 

+ 0.01 

- 0.01 

1.013 ± 0.027 

25 , 29 

0.032 


± 0.015 

+ 0.019 

- 0.004 

+ 0.003 

- 0.003 

0.888 ± 0.001 




0 < 

^ jetl ,2 ^ ^ 


14 , 17 

29.4 


± 0.8 

+ 0.4 

- 1.9 

+ 1.1 

- 1.1 

1.006 ± 0.013 

17 , 21 

12.7 


± 0.4 

+ 0.4 

- 1.4 

+ 0.5 

- 0.5 

1.006 ± 0.005 

21 , 25 

4.6 


± 0.3 

+ 0.2 

- 0.5 

+ 0.3 

- 0.2 

1.004 ± 0.014 

25 , 29 

1.2 


± 0.1 

+ 0.2 

- 0.1 

+ 0.10 

- 0.04 

0.972 ± 0.020 

29 , 35 

0.21 


± 0.04 

+ 0.07 

- 0.02 

+ 0.02 

- 0.01 

1.043 ± 0.038 


1 < 

< 2 . 4 , -1 < < 0 


14 , 17 

15.4 


± 0.4 

+ 2.3 

- 0.5 

+ 0.8 

- 0.5 

1.040 ± 0.019 

17 , 21 

6.0 


± 0.2 

+ 0.5 

- 0.1 

+ 0.3 

- 0.3 

1.039 ± 0.017 

21 , 25 

1.4 


± 0.1 

+ 0.2 

- 0.1 

+ 0.05 

- 0.04 

1.007 ± 0.019 

25 , 29 

0.38 


± 0.05 

+ 0.03 

- 0.18 

+ 0.03 

- 0.03 

1.031 ± 0.017 

29 , 35 

0.057 


± 0.015 

+ 0.013 

- 0.022 

+ 0.001 

- 0.001 

0.913 ± 0.013 



1 < 

< 2 . 4 , 0 < 

^ jet 2 ^ ^ 


14 , 17 

56.7 


± 0.8 

+1.6 

- 1.0 

+ 2.1 

- 2.1 

1.007 ± 0.024 

17 , 21 

26.6 


± 0.4 

+ 0.5 

- 0.5 

+ 1.3 

- 1.0 

1.021 ± 0.018 

21 , 25 

11.2 


± 0.3 

+ 0.1 

- 0.7 

+0.6 

- 0.5 

1.013 ± 0.014 

25 , 29 

3.7 


± 0.2 

+ 0.1 

- 0.2 

+ 0.1 

- 0.2 

0.999 ± 0.015 

29 , 35 

1.3 


± 0.1 

+ 0.05 

- 0.05 

+ 0.1 

- 0.1 

1.003 ± 0.004 

35 , 41 

0.35 


± 0.04 

+ 0.04 

- 0.04 

+ 0.03 

- 0.02 

1.026 ± 0.009 

41 , 48 

0.064 


± 0.016 

+ 0.007 

- 0.009 

+ 0.003 

- 0.005 

1.001 ± 0.004 

48 , 55 

0.019 


± 0.009 

+ 0.014 

- 0.005 

+ 0.001 

- 0.001 

1.000 ± 0.014 


Table 3: Measured cross section as a function of for events with < 0.75. The 
measurement is divided into six regions of the pseudorapidities of the jets. For further 
details, see the caption to Table [|. 
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bin ( GeV ) 

da/dE^^^^ Astat A ^ y^t Aes ( pb / GeV ) 

Ghad 

1 < < 2.4 


14 , 17 

78.4 

± 1.4 

+6.6 

- 5.4 

+ 2.3 

-2.8 

1.000 ± 0.024 

17 , 21 

38.2 

± 0.7 

+ 2.5 

- 0.9 

+ 1.7 

-1.6 

1.019 ± 0.025 

21 , 25 

16.9 

± 0.5 

+ 0.9 

- 0.3 

+0.8 

- 0.7 

1.016 ± 0.025 

25 , 29 

6.3 

± 0.3 

+ 0.4 

-0.1 

+ 0.3 

- 0.3 

1.012 ± 0.025 

29 , 35 

2.5 

± 0.1 

+0.1 

-0.0 

+0.1 

-0.1 

1.020 ± 0.032 

35 , 41 

0.92 

± 0.08 

+ 0.05 

-0.01 

+ 0.09 

- 0.06 

1.018 ± 0.031 

41 , 48 

0.31 

± 0.05 

+ 0.03 

- 0.04 

+0.02 

-0.02 

1.000 ± 0.026 

48 , 55 

0.12 

± 0.03 

+0.01 

-0.01 

+0.01 

-0.01 

0.995 ± 0.019 

55 , 65 

0.045 

± 0.016 

+ 0.014 

- 0.016 

+ 0.004 

- 0.003 

0.998 ± 0.003 


Table 3: (cont.) 
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^ jet2 


-^stat 

-^syst 

Aes (pb) 

Chad 



-1 < 

^jetl 

< 0 


-1.0, -0.5 

8.9 

±1.0 

+ 9.5 

- 1.4 

+ 0.5 

-0.8 

0.538 ± 0.021 

-0.5, 0.0 

88.2 

±2.8 

+ 15.6 

- 9.2 

+ 4.2 

- 3.2 

0.805 ± 0.031 

0.0, 0.5 

175.0 

±3.8 

+ 11.1 
- 15.7 

+ 5.8 

- 5.7 

0.881 ± 0.016 

0.5, 1.0 

195.8 

±4.0 

+ 6.7 

- 12.4 

+ 8.5 

- 6.3 

0.906 ± 0.012 

1.0, 1.5 

160.3 

±3.4 

+ 9.8 

- 5.0 

+ 6.5 

- 5.9 

0.921 ± 0.012 

1.5, 2.0 

130.0 

±3.3 

+ 6.9 

-2.6 

+ 4.2 

- 4.7 

0.913 ± 0.021 

2.0, 2.4 

89.3 

±3.5 

+ 5.1 

- 5.6 

+ 3.6 

- 3.1 

0.899 ± 0.021 



0 < 

^jetl ^ ^ 


-1.0, -0.5 

99.3 

±3.1 

+11.1 

- 1.9 

+ 4.7 

-2.8 

0.805 ± 0.012 

-0.5, 0.0 

272.0 

±4.6 

+ 10.7 
- 26.6 

+ 9.8 

- 9.1 

0.922 ± 0.010 

0.0, 0.5 

347.0 

±5.1 

+ 3.6 

- 37.8 

+12.1 

- 10.7 

0.952 ± 0.015 

0.5, 1.0 

347.1 

±5.0 

+2.2 

- 32.5 

+ 11.7 

- 9.2 

0.960 ± 0.013 

1.0, 1.5 

287.9 

±4.2 

+ 2.7 

- 19.5 

+ 9.3 

- 9.7 

0.971 ± 0.004 

1.5, 2.0 

208.1 

±3.7 

+1.6 

- 5.4 

+ 7.5 

- 7.2 

0.978 ± 0.012 

2.0, 2.4 

125.1 

±3.7 

+ 1.9 
- 7.9 

+ 3.9 

- 4.3 

0.967 ± 0.020 



1 < 7 ] 

jetl ^ 2.4 


-1.0, -0.5 

115.4 

±3.2 

+ 20.5 

-6.8 

+ 5.2 

- 4.8 

0.866 ± 0.022 

-0.5, 0.0 

241.7 

±4.2 

+6.6 

- 3.5 

+ 8.5 

- 8.3 

0.933 ± 0.011 

0.0, 0.5 

315.7 

±4.6 

+ 3.8 

- 13.3 

+ 11.4 

-11.2 

0.953 ± 0.018 

0.5, 1.0 

277.5 

±4.2 

+ 1.2 
- 18.6 

+ 8.5 

- 9.1 

0.992 ± 0.002 

1.0, 1.5 

104.0 

±2.4 

+ 4.1 

- 3.5 

+ 3.8 

- 4.0 

1.021 ± 0.007 

1.5, 2.0 

37.2 

±1.5 

+0.8 

- 1.4 

+ 1.5 

- 1.3 

1.018 ± 0.008 

2.0, 2.4 

13.2 

±1.1 

+ 0.9 

- 2.3 

+ 0.5 

- 0.4 

1.031 ± 0.017 


Table 4: Measured cross section as a function for events with > 0.75. 

The measurement is divided into three regions of the pseudorapidity of the other 
jet. For further details, see the caption to Table [J. 
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^ jet2 

da 

-^stat 

-^syst 

Aes (pb) 

Chad 



-1 < 

^jetl 

< 0 


-0.5, 0.0 

2.2 

±0.4 

+ 1.4 

-0.8 

+0.1 

-0.1 

1.088 ± 0.113 

0.0, 0.5 

15.8 

±1.0 

+ 2.9 

- 0.5 

+0.6 

-0.6 

1.003 ± 0.019 

0.5, 1.0 

40.0 

±1.7 

+ 2.7 

-2.2 

+ 1.3 

-1.0 

1.008 ± 0.008 

1.0, 1.5 

52.4 

±1.8 

+6.2 

- 0.4 

+2.6 

-1.8 

1.033 ± 0.027 

1.5, 2.0 

60.4 

±2.0 

+ 7.5 

- 4.6 

+ 3.0 

- 2.7 

1.028 ± 0.002 

2.0, 2.4 

54.1 

±2.4 

+ 6.9 

- 3.5 

+ 2.7 

- 1.4 

1.038 ± 0.024 



0 < 



-1.0, -0.5 

4.1 

±0.5 

+ 1.5 

-0.6 

+ 0.3 

-0.1 

1.046 ± 0.020 

-0.5, 0.0 

51.7 

±1.9 

+ 4.7 

-2.2 

+1.6 

-1.6 

1.004 ± 0.002 

0.0, 0.5 

131.7 

±3.1 

+ 3.3 

- 11.7 

+6.6 

- 5.2 

1.002 ± 0.005 

0.5, 1.0 

195.4 

±3.8 

+ 4.4 

- 14.4 

+ 7.4 

- 7.7 

1.002 ± 0.002 

1.0, 1.5 

230.1 

±3.7 

+ 4.8 

- 11.3 

+ 9.7 

- 9.4 

0.998 ± 0.023 

1.5, 2.0 

262.4 

±4.1 

+ 5.7 

- 4.1 

+ 11.5 

-10.2 

1.012 ± 0.014 

2.0, 2.4 

243.7 

±5.0 

+ 6.7 

- 7.1 

+10.1 

- 9.1 

1.020 ± 0.011 



1 <rj 

jetl ^ 2.4 


-1.0, -0.5 

20.5 

±1.2 

+ 5.1 

- 0.3 

+ 0.9 

-0.8 

1.056 ± 0.020 

-0.5, 0.0 

135.0 

±3.0 

+ 13.5 

- 5.1 

+ 6.9 

- 5.0 

1.028 ± 0.022 

0.0, 0.5 

288.8 

±4.4 

+ 5.6 

- 6.5 

+12.1 

-11.6 

1.016 ± 0.013 

0.5, 1.0 

403.6 

±5.2 

+ 9.6 

- 9.7 

+ 17.5 

- 15.5 

1.006 ± 0.018 

1.0, 1.5 

403.0 

±4.8 

+ 24.2 

- 15.6 

+ 16.3 

- 16.0 

1.008 ± 0.029 

1.5, 2.0 

355.2 

±4.5 

+ 27.6 

-20.8 

+ 13.6 

- 14.4 

1.013 ± 0.026 

2.0, 2.4 

305.0 

±5.3 

+ 27.4 

-11.1 

+11.0 

-10.8 

1.015 ± 0.021 


Table 5: Measured cross section as a function for events with < 0.75. 

The measurement is divided into three regions of the pseudorapidity of the other 
jet. For further details, see the caption to Table |^. 
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bin 

da / dxff ^ 

'^stat 

-^syst 

^ES ( pb ) 

Chad 

14 < < 17 GeV 

0.1, 0.2 

421.4 

±14.1 

+ 49.9 

- 33.5 

+ 19.5 

- 16.6 

0.970 ± 0.002 

0.2, 0.3 

631.6 

±16.2 

+ 48.2 

- 56.9 

+ 19.3 

- 25.7 

0.994 ± 0.027 

0.3, 0.4 

607.2 

±15.0 

+ 38.2 

- 14.3 

+20.2 

- 23.0 

0.992 ± 0.022 

0.4, 0.5 

599.9 

±14.4 

+ 36.7 

- 4.2 

+ 26.6 

- 19.9 

1.002 ± 0.023 

0.5, 0.6 

621.2 

±14.6 

+ 49.8 

- 6.5 

+21.6 

- 17.9 

1.016 ± 0.026 

0.6, 0.7 

695.8 

±15.4 

+ 33.2 

-12.0 

+ 24.6 

- 24.5 

1.035 ± 0.014 

0.7, 0.8 

842.1 

±16.5 

+ 97.5 

-6.6 

+ 33.4 

- 25.9 

1.113 ± 0.002 

0.8, 1.0 

1784.0 

±18.8 

+ 128.6 

- 78.1 

+ 57.5 

- 52.7 

0.906 ± 0.016 

17 < < 25 GeV 

0.1, 0.2 

251.2 

±9.6 

+11.2 

- 15.0 

+ 15.1 

- 13.3 

1.006 ± 0.019 

0.2, 0.3 

419.2 

±11.3 

+ 28.5 

- 20.3 

+ 20.3 

- 19.9 

1.005 ± 0.013 

0.3, 0.4 

523.1 

±12.5 

+ 39.1 

- 21.3 

+ 27.3 

-21.2 

1.008 ± 0.019 

0.4, 0.5 

539.7 

±12.5 

+ 10.3 

- 14.3 

+22.6 

- 23.2 

1.008 ± 0.012 

0.5, 0.6 

597.3 

±13.2 

+ 15.3 
- 18.1 

+ 28.0 

- 24.7 

1.025 ± 0.018 

0.6, 0.7 

641.3 

±13.6 

+ 21.3 

-12.6 

+ 29.4 

- 23.1 

1.022 ± 0.019 

0.7, 0.8 

812.2 

±15.2 

+21.2 

-21.0 

+ 34.0 

- 32.3 

1.079 ± 0.002 

0.8, 1.0 

1791.5 

±17.9 

+ 30.6 

- 154.1 

+ 65.9 

- 62.4 

0.933 ± 0.018 

25 < < 35 GeV 

0.1, 0.4 

40.0 

±1.8 

+ 8.4 

- 0.7 

+2.0 

-2.2 

1.002 ± 0.025 

0.4, 0.6 

89.8 

±3.3 

+ 5.4 

-2.1 

+ 3.9 

- 3.5 

1.011 ± 0.011 

0.6, 0.8 

130.0 

±4.0 

+ 2.3 

- 3.3 

+ 5.1 

- 5.7 

1.019 ± 0.021 

0.8, 1.0 

332.2 

±7.3 

+ 30.7 

- 32.0 

+11.6 

- 12.7 

0.948 ± 0.016 

35 < < 90 GeV 

0.1, 0.4 

3.9 

±0.6 

+ 0.4 

-2.0 

+ 0.3 

-0.2 

0.991 ± 0.018 

0.4, 0.6 

12.7 

±1.2 

+1.1 

- 0.3 

+1.1 

-0.8 

1.017 ± 0.020 

0.6, 0.8 

24.5 

±1.7 

+ 1.3 

- 2.5 

+2.2 

- 1.3 

1.017 ± 0.015 

0.8, 1.0 

77.6 

±3.6 

+1.6 

-6.1 

+ 4.5 

- 3.3 

0.961 ± 0.011 


Table 6: Measured cross section as a function of in four regions of 
For further details, see the caption to Table [^. 


25 
























Figure 1: Examples of (a) direct and (b) resolved dijet photoproduction diagrams 
in LO QCD. 
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Figure 2: Comparison of data and MC simulation for (a) (b) i/jb, (c) 

(d) , (e) and (f) The data are shown as points compared to Herwig 

(solid line) and Pythia (dashed line). Also shown is the LO component of direct 
photon processes in Herwig (hatched area). The simulated sample is normalised to 
the data and the fraction of direct and resolved photon processes combined according 
to a x^-fit lo the distribution in (a). 
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Figure 3: Measured cross sections as a function of |cos6**| for (a) < 0.75 

and (b) > 0.75, compared to NLO predictions. The data are shown with 

statistical errors (inner bars) and statistical and systematic uncertainties added in 
quadrature (outer bars). The uncertainty due to that of the jet energy scale is shown 
as the shaded band. The NLO prediction corrected for hadronisation effects is shown 
calculated using the GRV-HO and CTEQ5M1 PDFs for the photon and proton, 
respectively, and the scale set to Et/ 2 (solid line). The hatched band represents 
the quadratic sum of the theoretical uncertainties as discussed in Section [70.4 The 
prediction using the AFG-HO photon PDF is also shown (dashed line). In (c) 
the cross sections are area-normalised and the data shown for < 0.75 (solid 
points) and > 0.75 (open circles). 
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Figure 4: Measured cross section as a function of for events with > 0.75. 
The measurement is divided into six regions of the pseudorapidities of the jets. The 
cross sections are multiplied by the scale factor indicated in brackets so that all 
regions can be displayed in the same figure. For further details, see the caption to 

Fig. |. 
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Figure 5: Ratio of cross sections to the central theoretical prediction as a function 
of for events with x°’’® > 0.75. The measurement is divided into six regions 
of the pseudorapidities of the jets. For further details, see the caption to Fig. [^. 
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Figure 6: Measured cross section as a function of for events with < 0.75. 
The measurement is divided into six regions of the pseudorapidities of the jets. The 
cross sections are multiplied by the scale factor indicated in brackets so that all 
regions can be displayed in the same figures. For further details, see the caption to 

Fig. |. 
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Figure 7: Ratio of cross sections to the central theoretical prediction as a function 
of for events with < 0.75. The measurement is divided into six regions 
of the pseudorapidities of the jets. For further details, see the caption to Fig. [^. 
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Figure 8: Measured cross section as a function for events with > 0.75. 

The measurement is divided into three regions of the pseudorapidity of the other 
jet. For further details, see the caption to Fig. |^. 
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Figure 9: Measured cross section as a function for events with < 0.75. 

The measurement is divided into three regions of the pseudorapidity of the other 
jet. For further details, see the caption to Fig. 
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Figure 10: Measured cross section as a function of in four regions of 
compared to NLO predictions. For further details, see the caption to Fig. |^. 
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Figure 11: Ratio of cross sections to the NLO prediction using GRV-HO and 
CTEQ5M1 as the photon and proton PDFs, respectively, and the scale set to -E't/2 
as a function of in four regions of . The data are shown with statistical 
errors (inner bars) and statistical and systematic uncertainties added in guadrature 
(outer bars). The uncertainty due to that of the jet energy scale is shown as the 
shaded band. The theoretical uncertainty is shown as the hatched band. Predictions 
using MRST99 (dashed line) for the proton are also shown. 
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Figure 12: Ratio of cross sections to the NLO prediction using AFG-HO and 
CTEQ5M1 as the photon and proton PDFs, respectively, and the scale set to Et/^ 
as a function of in regions of . For further details, see the caption to 

Fig. [ 7 ]. 
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Figure 13: Measured cross section as a function of for a fixed range 

of transverse energy of the leading jet, 25 < < 35 GeV, compared to MC 


T 

simulation and NLO predictions for (a) 0 < < 1, (b) > 0.8 and (c) 

< 0.8. The typical magnitude of the statistical and systematic uncertainties 
added in guadrature is ±10%. The NLO prediction, corrected for hadronisation 
effects, calculated using the GRV-HO and GTEQ5M1 PDFs for the photon and 
proton, respectively, is shown as the thick solid line. The shaded band represents 
the guadratic sum of the systematic uncertainties, as discussed in Section \i0.^ . The 
prediction using the AFG-HO photon PDF is shown as the dashed line. The predic¬ 
tion o/Herwig, calculated using the GPV-LO and GTEQfL PDFs for the photon 
and proton, respectively, is normalised to the first point in the range 0 < < 1; 


it is shown as the histogram. 
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